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Abstract 

A dilute polymer solution is modeled as a suspension of non-interacting 
Hookean dumbbells and the effect of excluded volume is taken into account by 
incorporating a narrow Gaussian repulsive potential between the beads of each 
dumbbell. The narrow Gaussian potential is a means of regularising a delta- 
function potential — it tends to the delta-function potential in the limit of the 
width parameter /i going to zero. Exact predictions of viscometric functions 
in simple shear flow are obtained with the help of a retarded motion expansion 
and by Brownian dynamics simulations. It is shown that for relatively small 
non-zero values of fi, the presence of excluded volume causes a swelling of 
the dumbbell at equilibrium, and shear thinning in simple shear flow. On the 
other hand, a delta function excluded volume potential does not lead to either 
swelling or to shear thinning. Approximate viscometric functions, obtained 
by assuming that the bead-connector vector is described by a Gaussian non- 
equilibrium distribution function, are found to be accurate above a threshold 
value of n, for a given value of the strength of excluded volume interaction, z. 
A first order perturbation expansion reveals that the Gaussian approximation 
is exact to first order in z. The predictions of an alternative quadratic excluded 
volume potential suggested earlier by Fixman (J. Chem. Phys., 1966, 45, 785; 
793) are also compared with those of the narrow Gaussian potential. 



1 Introduction 



The fact that two parts of a polymer chain cannot occupy the same place at the 
same time due to their finite volume has been recognised in the polymer literature 
for many years now as being an extremely important microscopic phenomenon that 
governs the macroscopic behavior of polymer solutions. Like hydrodynamic 

interaction, the excluded volume effect influences the properties of polymer solutions 
even in the limit of extremely long chains because it is responsible for segments 
remote from each other along the polymer chain interacting with each other. 

While the effect of excluded volume on static properties of polymer solutions has 
been widely studied, there have been very few attempts at examining its influence 
on properties far from equilibrium. Excluded volume effects can be incorporated 
into bead-spring chain models for polymer solutions in a relatively straightforward 
manner by adding the excluded volume interaction force between a particular bead 
and all the other beads in the chain (pairwise) to the other potential forces that are 
acting on the bead. An noteworthy aspect of this approach is the kind of repulsive 
potential that is adopted to represent the excluded volume interactions. In static 
theories of polymer solutions, the excluded volume interaction is typically assumed 
to be a very short range 5-function potential. 

Fixmani and more recently Ahn et al.0 have attempted to predict the rheological 
properties of dilute polymer solutions by approximately incorporating the effects of 
both hydrodynamic interaction and excluded volume in a self-consistent manner 
into a bead-spring chain model. (Ahn et al. also include finitely extensible springs 
in place of Hookean springs). In order to obtain a solvable model, Fixmanu used 
a repulsive quadratic excluded volume potential in place of a ^-function potential. 
This leads to a tractable model since the bead-connector vectors are then described 
by a Gaussian non-equilibrium distribution function. Results obtained with the 
quadratic excluded volume potential have, however, not been compared so far with 
the results of other models for the excluded volume potential. 

Andrews et al.0 have recently carried out a numerical study of the influence 
of excluded volume interactions on rheological and rheooptical properties of dilute 
solutions, with the help of Brownian dynamics and configuration biased Monte Carlo 
simulations. A bead-spring chain model, with "Fraenkel" springs between beads and 
a Morse potential to represent excluded volume interactions, was used to to model 
the flexible polymer molecule. Attention was largely confined to the prediction of 
properties in elongational flow and transient shear flow. 

The predictions of their theories in the limit of long chains have not been consid- 
ered by Fixman,! Ahn et alE and Andrews et al.ll On the other hand, the universal 
character of excluded volume effects have been studied using renormalisation group 
theory methods based on kinetic theory models (with a 5-function excluded volume 
potential) by Ottinger and coworkers. 3il 

While the work of Andrews et al.@ is based on Brownian Dynamics simulations, 
the accuracy of the other approximate treatments of excluded volume cited above 
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has not been tested by comparison with Brownian Dynamics simulations (which are 
an ideal tool for testing approximations for nonlinear effects). This is in contrast to 
the situation that exists for kinetic theory models that only incorporate hydrody- 
namic interaction effects, where extensive comparisons between the exact results of 
Brownian Dynamics simulations and various approximations have been madeMi! 

It is the purpose of this paper to examine the influence of the excluded volume 
effect on the rheological properties of a dilute polymer solution by using a narrow 
Gaussian potential to describe the excluded volume interactions. Since the narrow 
Gaussian potential tends to the ^-function potential in the limit of a parameter fi 
(that describes the width of the potential) going to zero, it provides a means of 
evaluating results obtained with a singular 5-function potential. Compared to the 
^-function potential, analytical calculations are not significantly harder with the 
narrow Gaussian potential; quite often, upon setting /i = at the end of a cal- 
culation, the predictions of a 5-function potential can be obtained. Furthermore, 
since Brownian dynamics simulations cannot be performed with a 5-function poten- 
tial, simulations carried out with the narrow Gaussian potential for small values of 
the parameter [i provide a means of asymptotically obtaining the predictions of a 
5-function potential model. 

Any molecular theory that seeks to describe the dynamics of polymers in good 
solvents must simultaneously incorporate both the microscopic phenomena of hy- 
drodynamic interaction and excluded volume, since hydrodynamic interaction effects 
have been shown to have an unavoidable influence on the dynamic behavior of poly- 
mer solutions. However, it would be difficult to explore the consequences of such 
a theory for two reasons. Firstly, the incorporation of hydrodynamic interaction 
would lead to the complication of multiplicative noise. Secondly, since Brownian 
dynamics simulations for long chains would be extremely computationally intensive, 
any approximations that are developed can only be tested for very short chains. For 
these reasons, and being in the nature of a preliminary investigation, we examine 
excluded volume effects independently from hydrodynamic interaction effects, and 
confine attention here to a Hookean dumbbell model for the polymer. This enables 
the careful evaluation of various approximations. It is hoped that, in the future, 
the best approximation can be used for chains with both hydrodynamic interaction 
and excluded volume, with the ultimate objective of calculating universal properties 
predicted in the long chain limit. 

It must be noted that unlike in the case of the quadratic potential, the use of the 
narrow Gaussian potential does not lead to exact solvability. Indeed, as has been 
observed in earlier treatments of other non-linear microscopic phenomenon such as 
hydrodynamic interaction and internal viscosity, it is found that it is not possible to 
obtain an analytical solution valid at all shear rates. As a result, retarded motion 
expansions, Brownian dynamics simulations, and perturbative and non-perturbative 
approximation procedures are used in this work to obtain the material functions 
predicted by the narrow Gaussian potential. 

An important consequence of using the narrow Gaussian potential is that the 
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nature of Fixman's quadratic excluded volume potential can be explored. Predic- 
tions of Fixman's quadratic excluded volume potential in simple shear flow will be 
compared with the predictions of the narrow Gaussian potential. 

The outline of this paper is as follows. In section 2, the basic equations re- 
quired to discuss the dynamics of Hookean dumbbells in good solvents are derived. 
In section 3, the various material functions considered in this paper are defined. 
A retarded motion expansion for the stress tensor, for arbitrary excluded volume 
potentials, is derived in section 4, and power series expansions for the material 
functions in simple shear flow are obtained. Section 5 is devoted to examining the 
consequences of describing the excluded volume interaction with a narrow Gaussian 
potential. In section 5.1, expressions for the material functions predicted with this 
potential at zero shear rate are obtained by using the retarded motion expansion. 
The formulation of a Brownian dynamics simulation algorithm is discussed in sec- 
tion 5.2, and a Gaussian approximation for the configurational distribution function 
is introduced in section 5.3. In section 5.4, a first order perturbation expansion 
in the strength of the excluded volume interaction is derived. Fixman's theory for 
dumbbells is presented in section 6, in terms of the framework adopted for the rest 
of the discussion in this paper. The results of the various exact and approximate 
treatments are compared and discussed in section 7, and the principal conclusions 
of the paper are summarised in section 8. 

2 Basic Equations 

The Hookean dumbbell model represents a macromolecule by a mechanical model 
that consists of two identical beads connected by a spring. The solvent in which 
the beads are suspended is assumed to be Newtonian, and attention is restricted to 
flows which have a homogeneous velocity field, ie. of the form v = v + K,(t) ■ r, 
where Vq is a constant vector, n(t) is a traceless tensor, and r is the position vector 
with respect to a laboratory-fixed frame of reference. The instantaneous position of 
the beads are specified by bead position vectors r\ and r 2 . 

This paper examines the consequence of introducing an excluded volume inter- 
action between the beads of the dumbbell, so that the total potential experienced by 
the beads, <j), is the sum of the spring potential S, and the excluded volume potential 
E. The force on bead v due to this potential, F^f \ is then given by F^f 1 = — d<p/ dr u . 
For Hookean springs, the spring potential is given by S = (1/2) HQ 2 , where H is 
the spring constant. Two forms of the excluded volume potential are considered in 
this work, namely, the narrow Gaussian potential, and Fixman's quadratic poten- 
tial. These are discussed in greater detail subsequently. Here, we summarise the 
governing equations that are valid for an arbitrary choice of the excluded volume 
potential E. 

For homogeneous flows, the configurational distribution function tp(Q, t) depends 
only on the internal configuration of the dumbbell, specified by the bead-connector 
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vector Q = r 2 — r 1; and not on the center of mass. The quantity i[)(Q,t) dQ is 
then the probability that at time t the dumbbell has a configuration in the range 
Q to Q + dQ. Using the framework of polymer kinetic theory! one can show that 
the distribution function ip(Q,t), in the presence of excluded volume, satisfies the 
following diffusion equation, 

d^P d J 2 d<f> 2k B T 

where, ( is the bead friction coefficient (ie., ( = 67n] s a for spherical beads with 
radius a, in a solvent with viscosity r] s ), k B is Boltzmann's constant, and T is the 
absolute temperature. 

The stress tensor, r, in a polymer solution is considered to be given by the 
sum of two contributions, r = t s + t p , where t s is the contribution from the 
solvent, and t p is the polymer contribution. Since the solvent is assumed to be 
Newtonian, r s = — n s 7, where 7 is the rate of strain tensor, 7 = (Vi>)(t) + (Vi^(t). 
The rheological properties of a dilute polymer solution may thus be obtained by 
calculating the polymer contribution to the stress tensor, t p . For a dumbbell model 
in the presence of excluded volume, it is given by the Kramer's expression,! 

r p = -n(Q^) + nk B Tl (2) 

Here, n is the number density of polymers, and angular brackets represent averaging 
with respect to the configurational distribution function ip(Q,t). 

For both the excluded volume potentials considered here, it will turn out that 
calculation of the second moment (QQ) is necessary in order to evaluate the average 
in equation A time evolution equation for the second moment can be obtained by 
multiplying the diffusion equation (|IJ) by QQ and integrating over all configurations, 



j f (QQ) = k ■ (QQ) + (QQ) .^ + ^1-1 



(3) 



It proves convenient for subsequent calculations to introduce the following di- 
mensionless variables, 



t 




Q* = J—Q, k* = \ h k, <f>* = S* + E* (4) 



where, Xh = (C/^H) is the familiar time constant, S* = S/ksT = (1/2) Q* 2 
and E* = Ejk^T are the non-dimensional Hookean spring potential and the non- 
dimensional excluded volume potential, respectively. Note that ^-solvent values of 



a typical time scale (A#) and a typical length scale (\Jk B T / H) are used for the 
purpose of non-dimensionalisation, regardless of the form of the excluded volume 
potential. 
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The intramolecular force is expected to be in the direction of the bead-connector 
vector. Therefore, it can be written in terms of non-dimensional variables as, 

H*(Q*)Q* (5) 



dQ* 

where H*(Q*) is an arbitrary function of the magnitude of Q*. As a result, the 
diffusion equation ([I]), in terms of non-dimensional variables, is given by, 

°* a r K ..Q._V (Q . )cr U + i » (6) 



dt* dQ* I ^ 2 v ^ y ^ J r 2 <9Q* <9Q 
while, the Kramers expression (||) assumes the form, 

T ' -(H*(Q*)Q*Q*) + 1 (7) 



nksT 

and the second moment equation (^|) becomes, 

A(Q*Q*) = K * . (Q*Q*) + (Q*Q*) • K * T - (H*(Q*) Q*Q*) + 1 (8) 

Equation (||) is in general not a closed equation for the second moments since it 
depends on the form of the function H*(Q*) (ie. on the choice of excluded volume 
potential). 

On examining equations (j?]) and it is straightforward to see that the second 
moment equation (§) is nothing but the Giesekus expression for the stress tensor, 

= ^(Q*Q1 - *• ■ (Q*Q*> - (Q*Q*> ■ ^ T (9) 

While the equations derived in this section are valid for arbitrary homogeneous 
flows, in this paper we confine attention to the prediction of rheological properties 
in simple shear flows, defined in the section below. 



3 Simple Shear Flows 
3.1 Steady simple shear flow 

Steady simple shear flows are described by a tensor k which has the following matrix 
representation in the laboratory-fixed coordinate system, 

/0 1 0\ 

k = 7 (10) 
\0 0/ 

where 7 is the constant shear rate. 

The three independent material functions used to characterize such flows are 
the viscosity, r) p , and the first and second normal stress difference coefficients, 
^1 and ^2, respectively. These functions are defined by the following relations,^ 

T P xy = ~iV P ; r p xx - rl y = - 7 2 ; r p „ - 7* = -f ^2 (11) 
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3.2 Small amplitude oscillatory shear flow 



A transient experiment that is used often to characterise polymer solutions is small 
amplitude oscillatory shear flow, where the tensor n(t) is given by, 

/0 1 0\ 

K (t) = 70 cos cut (12) 

\o o oj 

Here, 70 is the amplitude, and u is the frequency of oscillations in the plane of 
flow. The polymer contribution to the shear stress, t*} x , depends on time through 
the relation,^ 

T yx = ~ 7 l'( UJ ) 7o cos uot - rj"(u) 70 sin tut (13) 

where 7/ and rj" are the material functions characterising oscillatory shear flow. They 
are represented in a combined form as the complex viscosity, rj* = rj' — i 77". 

In the linear viscoelastic flow regime, the stress tensor is described by the linear 
constitutive relation, 



(t) = - I dsG(t-s)i(s) (14) 

J — 00 

where G(t) is the relaxation modulus. As a result, for oscillatory shear flows with 
a small amplitude 70, expressions for the real and imaginary parts of the complex 
viscosity can be found in terms of the relaxation modulus from the expression, 

rf = / G(s)e- iujs ds (15) 
J 

Note that the zero shear rate viscosity r] Pi o arid the zero shear rate first normal 
stress difference ^i.o, which are linear viscoelastic properties, can be obtained from 
the complex viscosity in the limit of vanishing frequency, 

r/ p ,o = lim v '(u) ; * ll0 = Hm ^M- (16) 



4 Retarded Motion Expansion 

A retarded motion expansion for the stress tensor, derived previously for the FENE 
dumbbell model,! can be adapted to the present instance by recognising that the 
FENE spring potential (as indeed any choice of excluded volume potential) is but a 
particular example of a connector force potential between the beads of the dumbbell. 
In this section, we briefly summarise the development of a retarded motion expansion 
for an arbitrary choice of the excluded volume potential. Details of the derivation 
may be found in Bird et al.i 

We seek a solution of the diffusion equation (Rj) whereby the configurational 
distribution function, ip(Q,t), can be written as a product of an equilibrium contri- 
bution and a flow contribution, 

V(Q,*)=V'eq(Q)0a(O,t) (17) 
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The governing equation for the flow contribution fl (Q,t), 



dt* XdQ* dQ* } ' K 2 I dQ* ' OQ* OQ* ' OQ* J 1 ' 

can be obtained by substituting equation (|I7|) into the diffusion equation ([]), and 
exploiting the fact that the equilibrium distribution function is given by, 

^e q (Q) =AT eq e-^ (19) 

where A/" eq is the normalisation constant. 

Regardless of the form of the excluded volume potential, at steady state, an 
exact solution to equation flT8| ) can be found for all homogeneous 'potential flows.! 
For more general homogeneous flows, however, it is necessary to seek a perturbative 
solution. The flow contribution, 0a (Q,t), is assumed to be expandable in a power 
series in the velocity gradients, 

<MQ,t) = 1 + 01 + 02 + 03 + --- (20) 

where <frk is of order k in the velocity gradient. 

Partial differential equations governing each of the <fik may be obtained by substi- 



tuting equation ( P0|) into equation (|1 
the procedure suggested in Bird et al., 



and equating terms of like order. Following 
one can judiciously guess the specific forms 
for the functions (f>k by noting that each of these functions must have certain prop- 
erties. The form of the function (f>i which satisfies these requirements is, 



2 

while the form of 02 can be guessed to be 



i = \Q*-i-Q* (21) 



^2 



(Q* • 7 ■ Q*f ~ ^ ( Q* 4 ) eQ t r (7 • 7) + A *(Q*) Q* 7 " Q* (22) 



60 



' cq 



where, ( ) denotes an average with ip e q(Q), <*> is the vorticity tensor, defined here 
as a; = k* — k* t , and the scalar function A*(Q*) obeys the following second order 
differential equation, 

+ ( ^ - H*(Q*) Q*)^-2 A*(Q*) H*(Q*) = 1 (23) 

It is difficult to suggest boundary conditions for equation (|23| ) other than to say that 
the solution must be such that tp(Q,t) is bounded. In the case of FENE springs, 
it is possible to explicitly obtain the particular solution of a similar second order 
differential equation. 

It is clear from equation (Q) that at steady state, the stress tensor can be found 
once (Q*Q*) is known. The second moment (Q*Q*) can be found correct to second 
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order in velocity gradients by using the power series expansion ( ^0[) for 0a (Q, t), and 
the specific forms for (f>i and $2 in equations fl2"T| ) and (^), respectively. This leads 
to the following expression for the stress tensor, correct to third order in velocity 
gradients, 



A " ' H '(Q 2 ) cq {7} + #f^ 2 (Q 4 ) cq {27 2 -(7---^-7)} 



nk B T 3 \k B T J /cqL,J 30 \k B Tj /eq 

+ ^|(^) 2 W 4 ^)e q {(7 2 ^-^-7 2 ) + ^-(7-^-^-7) 

- (7 • a> - w • 7) • a> } + . . . (24) 

The Cayley-Hamilton theorem has been used to eliminate the term 7 3 in equa- 
tion (p4|) , and an isotropic term that does not affect the rheological properties has 
been dropped. 

The stress tensor in simple shear flow, for small values of the non-dimensional 
shear rate Xnj, can be found by substituting equation ( |10"D for the rate of strain 
tensor, in equation Using the definitions of the viscometric functions in equa- 
tion (pi]) , the following power series expansions are obtained, 



eq 



\ H nk B T 3 \k B Tj ^ /oq I 15 \k B T J ^ /eq 70 \k B T 
1 / H 



90 \k B T 
*i 2 / H ^ 2 



) (Q%(Q 2 ) cq }(^i) 2 + ... (25) 



\ 2 H nk B T 15 U B 77 ^ /eq 



(Q 4 )eo + • • • (26) 



Clearly, equation ( pq ) indicates that one must expand to higher orders in velocity 
gradients before the shear rate dependence of the first normal stress difference can be 
obtained. Zero shear rate properties, however, can be obtained from equations (|2"5| ) 
and (EH). 



5 The Narrow Gaussian Potential 

In the static theory of polymer solutions it is common to represent the dimensionless 
excluded volume potential, between two points on the polymer chain separated by 
a non-dimensional distance Q*, with the Dirac delta function, 

E* (Q*) = (2tt) 3/2 2 5(Q*) (27) 
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where, z = v ( H/27rksT ) 3 / 2 is a non-dimensional parameter which represents the 
strength of the excluded volume interaction, and in which v — which has the dimen- 
sions of volume — is called the 'excluded volume parameter'.! The parameter z is 
frequently used in theories that incorporate excluded volume, as it is considered to 
be the appropriate parameter to be used in perturbation expansions. As mentioned 
earlier, excluded volume interactions are taken into account in this work by means 
of a narrow Gaussian potential.^ The narrow Gaussian potential has the following 
form in terms of non-dimensional variables, 

= (-5?) (28) 

It is clear from equation (|28| ) that the non-dimensional parameter /i controls 
the extent of the excluded volume interaction, and as /i — > 0, the narrow Gaussian 
potential tends to the 5-potential. The narrow Gaussian potential, as mentioned 
earlier, serves as a means of regularising the singular 5-potential and consequently, 
permits the evaluation of results obtained with a (^-potential. 

With excluded volume interactions described by the narrow Gaussian potential, 
the function H*(Q*), which appears in equation (||]) for the non-dimensional force 
between the beads of the dumbbell, is given by, 

H*(Q*) = Hq(Q*) = 1 - (JL\ exp 

The complex form of this function implies that the diffusion equation (^) cannot 
be solved exactly analytically to obtain the non-equilibrium configurational distri- 
bution function ip(Q,t). Furthermore, the time evolution equation for the second 
moments (|8]) is not a closed equation for the second moments since it involves the 
higher order moment (Hq(Q*) Q*Q*) on the right hand side. As a result, pertur- 
bative methods, non-perturbative approximation procedures or numerical schemes 
must be used in order to obtain the material functions predicted by the narrow 
Gaussian potential. 

5.1 Zero shear rate properties 

The viscosity and the first normal stress difference predicted by the narrow Gaussian 
potential at low shear rates can be obtained from the equations (|25|) and (^6|), respec- 
tively, once the equilibrium averages that occur in these expressions are evaluated. 
For the narrow Gaussian potential, the equilibrium distribution function is given by 
equation (|1|), with the excluded volume contribution to the non-dimensional po- 
tential 0* given by equation (p8|). We denote the various non-dimensional moments 
of the equilibrium distribution for a narrow Gaussian potential by, 

1m = (j^) (Q 2m ) eq ; m = 1,2,3,... (30) 
10 



Q 



:2 



2^ 



(29) 



In order to obtain the viscosity at non-zero shear rates, it is necessary to find the 
function A* that satisfies the second order differential equation [with H*(Q*) = 
Hq(Q*)]. Unlike in the case of an FENE dumbbell, it has not been possible to 
obtain the particular solution to equation (P3[). As a result, attention is confined 
here to obtaining the zero shear rate predictions of the narrow Gaussian potential, 

Vp,o _ 1 



Xunk^T 3 
*i,o S 



qi (31) 

Q2 (32) 



X 2 H nk B T 15 

for which only the moments q\ and q2 are required. Alternative methods will be 
used in sections 5.2 to 5.4 to obtain the shear rate dependence of the viscometric 
functions. 

The non-dimensional equilibrium moments q\ and q<i can be obtained exactly, as 
will be shown in section 5.1.1 below. The need to calculate equilibrium moments is 
also frequently encountered in static theories of polymer solutions. In these theories, 
as mentioned earlier, it is common to represent excluded volume interactions with 
a delta-function excluded volume potential, and furthermore, to obtain universal 
predictions by considering the limit of long chains. The Hamiltonian then typically 
has two singular objects — making it impossible to evaluate equilibrium moments 
exactly. The most succesful approach so far towards approximately evaluating these 
moments has been to develop a perturbation expansion in the parameter z, and to 
use renormalisation group methods to refine the results of the perturbation calcula- 
tion.Bi An alternative and simpler non-perturbative route is the uniform expansion 
model.i The use of the narrow Gaussian potential — albeit in the simple context 
of Hookean dumbbells — provides an opportunity to compare the results of these 
approximate models with the exact solution. The rigorous solution is discussed in 
section 5.1.1 below, while the perturbation expansion and the uniform expansion 
models are discussed in sections 5.1.2 and 5.1.3, respectively. 



5.1.1 Exact solution 

It is straight forward to show that the moments q m are given by the ratio of two 
integrals, q m = (I m /I ), where, 

POO 

Ij= Q* 2l+2 exp {- (1/2) Q* 2 -E*}dQ*; j = 0,1,2,... (33) 

J 

In the limit of // — > and ji — > oo, these integrals can be evaluated analytically. 
Consider the quantities, 

which are the integrands for the integrals Ij, and 

^•(g*)^g^ +2 ex P {-ig* 2 } 
n 



Now, Pj(0) = Rj{0) = 0, for all values of p. At any non-zero value of Q*, it is clear 
from equation Q2"8|) that, 



Pj(Q*) ~^ Rj(Q*) as a* — * o or 00 

In other words, for all values of Q*, the quantities Pj{Q*) tend pointwise to 
Rj(Q*) as // tends to zero or to infinity. Furthermore, for all values of p, it can be 
shown that Pj(Q*) are bounded functions of Q* on [0,oo[. It then follows from a 
theorem of the calculus! that, 

C Q* 2m+2 exp{-iQ* 2 }rfQ* 

Jo Q exp{-ig* }c/g* 

As a result, the asymptotic values of q m for p — > and /i — > 00 are found to be 
independent of 2, and are equal to the ^-solvent values, 

Qi = 3 ; <? 2 = 15 ; g 3 = 105 ; . . . 



This implies — from equations ([H|) and ([32]) — that the use of a delta-function poten- 
tial to represent excluded volume interactions leads to the prediction of zero shear 
rate properties in good solvents which are identical to those in ^-solvents. 

Away from these limiting values of p, i.e. at non-zero finite values of p, the 
integrals Ij can be found by numerical quadrature. Here they have been evaluated 
using a routine given in Numerical recipesS for the integration of an exponentially 
decaying integrand. Discussion of zero shear rate property predictions in this case 
is taken up in section 7. 



5.1.2 Perturbation expansion 

Static theories for polymer solutions indicate that accounting for excluded volume 
interactions with a delta function potential leads to the prediction of a swelling (i.e. 
an increase in the mean square end-to-end distance) of the polymer chain, which 
is in close agreement with experimental observations.^'! In the case of Hookean 
dumbbell, however, we have seen above that the use of a delta function potential 
to account for the presence of excluded volume (which corresponds to the limit 
p —* 0), does not lead to any change in the prediction of equilibrium moments when 
compared to the ^-solvent case. It is worthwhile therefore to examine the nature 
of the perturbation expansion in z, and to compare it with the results of the exact 
calculation. 

Upon expanding e~ E * in a power series, the integral Ij has the form, 

/>oo 00 

/, = / dQ* u n(Q*) ; J = 0,1,2,... (34) 

where, 

u n {Q*) = K —r- Q* 2J+2 e" 5 Q * 2 ( E*(Q*) ) n (35) 
12 



In order to carry out a term by term integration of the functional series J2^=o u n(Q*) 
in equation Q3"3j), it is necessary for the series to be uniformly convergent on [0, oo[. 
For all values of z, and fi 7^ 0, uniform convergence can be established with the help 
of the Weierstrass M test.! Therefore, a term by term integration in equation fl3*3| ) 
can be caried out, except when /i = 0. Assuming that /j 7^ 0, and performing the 
integration in equation fl3"4|), one obtains, 



h = 2^ T(j + |) £ - ; 3 = 0,1,2,... (36) 

Consider the moment q\, which is undoubtedly the most interesting physical 
moment. Using the perturbation expansion for Ij, qi is given by the ratio, q\ = 
3 (S\/ So), where, So and Si are defined by, 

■l) n f z\ n /i 3 



° to n\ ( n + ^) 

(-l) n /z\ n /i 5 

to ™ ! v^J (^T^) 



* = E ^ ^ (37) 



If only the first order perturbation term is retained, one obtains, 

* = 3 ( 1 + (TW) (38) 

Curiously, the limit \i — > can be carried out in this case. It leads to a result 
which is in line with static theories, which indicate a finite non- vanishing effect due 
to the presence of (^-potential excluded volume interactions. However, such a limit 
cannot be carried out if higher order terms are retained, since both the sums So 
and Si diverge as fi — > 0. In static theories, higher order perturbation expansions 
are obtained by dropping divergent terms, as they are postulated not to matter. 
These divergent terms arise from products of the (^-potential for the same pair of 
interacting beads. In dumbbells, these are the only kind of beads present. As a 
result, there is no meaningful or mathematically consistent way of going beyond 
first order perturbation theory for dumbells in the limit \x — ► 0. 

Note that both the series, So and Si, are alternating series. Furthermore, for 
given values of z and [i, the terms decrease monotonically for large enough values of 
n. It follows then, from the Leibnitz criterion for alternating series,! that both So 
and Si converge for all values of z and /x 7^ 0. This suggests that, even though it is 
not possible to switch the integral and summation in equation (|34]) for fi = 0, the 
value of qi at n = can be found by setting it equal to the limit of qi(fi) as \i — » 0. 
We shall see in section 7 that such a limiting process is infeasible since it becomes 
numerically impossible to evaluate the sums ([37]) for small enough values of \i. 
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5.1.3 Uniform expansion model 

The uniform expansion model seeks to approximate the average (X) of any quan- 
tity X(Q), with (X) , where ( ) denotes an average with the Gaussian equilib- 
rium distribution function, 

<(Q)=A4' q exp{- Ag 2 } (39) 

with A/" c ' q = [3/27r6' 2 ] 3 / 2 . The aim is to find the parameter b' that leads to the 
best possible approximation. As may be expected, this depends on the quantity 
X(Q) that is averaged. The motivation behind the uniform expansion model, and 
details regarding the calculation of b' are given in appendix A. Since the equilibrium 
distribution function in the absence of excluded volume is Gaussian, this assumption 
expects the equilibrium distribution to remain Gaussian upon the incorporation of 
excluded volume, albeit with an increased end-to-end vector. 

If we define u m , such that u m = {H/k-QT) m (Q 2m )' eq , then clearly u m is the 
uniform expansion model approximation for q m . The uniform expansion model 
predictions of the zero shear rate properties are given by (|3T|) and ([32]), with q m 
replaced by u m . Results of material properties obtained by numerical quadrature, 
and by the uniform expansion model, are discussed in section 7. 



5.2 Brownian dynamics simulations 

Development of the retarded motion expansion has proved useful in obtaining exact 
expressions for the zero shear rate properties. At non-zero shear rates, exact pre- 
dictions of the viscometric functions can be obtained by solving the Ito stochastic 
differential equation, 



dQ* 



Q* dt + dW (40) 



which corresponds to the non-dimensional diffusion equation (||), and in which W 
is a three-dimensional Wiener process. 

For the narrow Gaussian potential, since H*(Q*) = Hq(Q*), equation fl4"0| ) is 
non- linear. As a result, it cannot be solved analytically. Two different Brownian 
dynamics simulation algorithms have been adopted here for the numerical solution 
of equation (|40|). Both schemes use a second order predictor-corrector algorithm 
with time-step extrapolation. El The first scheme obtains steady-state expectations 
by the simulation of a single long trajectory, and is based on the assumption of 
ergodicity.0 It has been used to obtain results at equilibrium, and for large val- 
ues of the shear rate. A second algorithm — which employs a variance reduction 
procedure — has been used at low values of the shear rate, since the variance for the 
viscometric functions is found to be relatively large at these shear rates. Reduc- 
tion in the variance is obtained by following a scheme suggested by Wagner and 
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Ottinger.Ea The scheme — which constructs an ensemble of trajectories over several 
relaxation times, from start-up of flow to steady-state — essentially consists of sub- 
tracting the rheological properties obtained from a parallel equilibrium simulation. 
While this doesn't change the average values of the properties, it significantly re- 
duces the fluctuations, since the fluctuations are virtually the same at zero and small 
shear rates. The results of these simulation algorithms are discussed in section 7. 



5.3 The Gaussian approximation 

The main obstacle (in the configuration space of the dumbbell) to obtaining the 
rheological properties predicted by a narrow Gaussian potential is that the second 
moment equation is not a closed equation. A closure problem has also been en- 
countered earlier in treatments of the phenomenon of hydrodynamic interaction and 
internal viscosity, where it has been shown that an accurate approximation can be 
obtained by assuming that the non-eguilibrium configurational distribution function 
is a Gaussian distributionSBIll'Bil In this section, a similar systematic approxi- 
mation procedure for the treatment of excluded volume interactions described by a 
narrow Gaussian potential is introduced. 

The assumption that if)(Q,t) is a Gaussian distribution, 

*) = Trwi i * ex P \-\q- (QQ)' 1 ■ Q 1 ( 41 ) 

( 2 *T /2 v /det(QQ) 1 2 > 

makes the second moment equation (^j) a closed equation, since the higher order 
moment (Hq(Q*) Q*Q*) can be expressed in terms of the second moment. On 
performing this reduction, it can be shown that the Gaussian approximation leads 
to the following closed second moment equation, 

<! 'Q*Q*) = n* ■ (Q*Q*) + (Q*Q*) ■ k* t - (Q*Q*) + I II + 1 

' det [ (Q*Q*) + ii 2 1] 

(42) 



dt 
where, 



n= [(Q*Q*) + fi 2 l]~ 1 ■ {Q*Q*} 



It is also straight forward to show that on introducing the Gaussian approxima- 
tion, the Giesekus expression for the stress tensor has the form, 

-(Q*Q*) + . = n+ 1 (43) 



nk sT ' 1 Jdet[(Q*Q*)+ n 2 l 



The steady state viscometric functions [defined by equations flTT|) 1 can therefore be 
found once equation fl42|) is solved for (Q*Q*). 

It is worth examining the nature of the polymer contribution to the stress tensor. 
In the limit // — > 0, II reduces to 1, and as a result, the presence of excluded volume 
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only makes an indirect contribution through its influence on the second moment 
(Q*Q*)- This follows from the fact that an isotropic contribution to the stress 
tensor makes no difference to the rheological properties of the polymer solution. On 
the other hand, for non-zero values of /i, the rheological properties are also affected 
directly by excluded volume. 

Under the Gausssian approximation, linear viscoelastic properties can be ob- 
tained by deriving a first order codeformational memory integral expansion for 
the stress tensor. The tensor (Q*Q*) is expanded, in terms of deviations from 
its isotropic equilibrium solution, upto first order in velocity gradient, 

(Q*Q*) =a 2 {l + e + ...) (44) 

where, the parameter a (commonly called the swelling ratio) is defined by, 

or = . . (45) 



Q 2 k 



eq 



Here, (Q 2 ) = (3k B T/ H), is the mean square end-to-end distance in the absence 
of excluded volume. Clearly, a represents the equilibrium swelling of the polymer 
chain caused by the presence of excluded volume, a is not an independent param- 
eter since at equilibrium the second moment equation ( f4"2| ) reduces to the following 
consistency relationship between z, /i and a, 



, 5 



z = (1 - a ~ z ) [a 2 + /^]3 (46) 

The well known scaling relation for the end-to-end distance with the number of 
monomers N, namely, J~{ Q 2 ) cq ~ iV 3//5 , may be obtained from equations fl45|) and 
( [46"D in the limit of large iV by noting that — since excluded volume is a pairwise 
interaction — z must scale as \JN for dumbbells.0 

Substituting the expansion ( [4*4] ) into the evolution equation ( f£2D leads to, 

4- e = k* + k* T - — e (47) 
dt* t* v ' 



where, 



z fi 2 



(46 

a 2 + fi 2 ) 7 / 2 . 

Furthermore, the stress tensor (|4*3D upto first order in the velocity gradient (without 
the Theologically unimportant isotropic contribution) is given by, 

4^ = -« e ( 49 ) 

where, 

H = a 2 {T*)- 1 (50) 
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Upon integrating equation (f47|) , which is a first order ordinary differential equa- 
tion for e, and substituting the result into equation (P?D, the following codeforma- 
tional integral expansion is obtained, 

T P(t) = -f t dsnk B TG{t-s)j {1] {t,s) (51) 

where ■ym is the codeformational rate-of-strain tensor! and the memory function 
G(t) is given by 

G(t) = He- (1/AflT,) (52) 

The product A#t* is usually interpreted as a relaxation time, and 7i as a relaxation 
weight. Clearly, the incorporation of excluded volume effects increases the relaxation 
time in a good solvent relative to a theta solvent by a factor of r*. 

The memory function G(t) can now be used to derive the linear viscoelastic 
material properties. Substituting G(t) = nk B TG(t), into equation (|15|), leads to, 



t*H t]"(uj) ujt* 2 H 



\ H nk B T 1 + {\ h t*uj) 2 ' \ 2 H nk B T 1 + (\ h t* uj) 2 
Upon taking the limit of uj — > we obtain, 

Vp,0 _ *nj . ^1,0 r,_*2. 



(53) 



, r*H ; ^ =2t**H (54) 

At moderate to large values of the shear rate, it is not possible to obtain ana- 
lytical expressions for the shear rate dependence of the viscometric functions, and 
consequently a numerical procedure is required. Since the second moment (Q*Q*) 
shares the symmetry of the flow field in simple shear flow, its Cartesian components 
can be denoted by, 

/ Si S4 \ 



(Q*Q*) 



s 4 s 2 
\0 s 3 J 



(55) 



Upon substituting equation ( p5|) into equation (|42|), a system of four first order or- 
dinary differential equations for the quantities Sj, j — 1, ... ,4 is obtained. Steady 
state viscometric functions, as functions of shear rate, can then be found by numeri- 
cally integrating these equations with respect to time (using a simple Euler scheme) 
until steady state is reached. Results obtained by this procedure are discussed in 
section 7. 



5.4 First order perturbation expansion in z 

The influence of excluded volume effects on the universal shear rate dependence of 
viscometric functions has been studied, as mentioned earlier, by using renormalisa- 
tion group methods. The renormalisation group theory approach is essentially 
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a method for refining the results of a low order perturbation expansion in z, by in- 
troducing higher order interactions effects so as to remove the ambiguous definition 
of the bead size. Results of this approach, based on a 5-function excluded volume 
potential, indicate that the presence of excluded volume has a non-trivial influence 
on the shear rate dependence of the viscometric functions. We have seen earlier 
in this work that as far as equlibrium swelling and zero shear rate properties are 
concerned, the use of a 5-function excluded volume potential leads to trivial results. 
In this section, a first order perturbation expansion in z — with a narrow Gaussian 
excluded volume potential — is constructed, in order to compare its predictions of 
shear rate dependence with those obtained with Brownian dynamics simulations, 
and with the Gaussian approximation. The dependence of the predictions on the 
width parameter \x is of particular interest. 

A first order perturbation expansion can be constructed, following the procedure 
suggested in references. EE from the second moment equation (||). We assume that 
the configurational distribution function -0 can be written as ipe + ip z , where ipg is the 
distribution function in the absence of excluded volume, i.e. in a ^-solvent, and ip z is 
the correction to first order in the strength of the excluded volume interaction. The 
averages performed with these contributions will be denoted by ( • • • ) g and ( • • • ) z , 
respectively. 

On equating terms of equal order, equation (||) can be rewritten as two equations, 
namely, a zeroth order second moment equation and a first order second moment 
equation. The zeroth order equation — which is linear in the moment (Q*Q*) g — is 
the well known second moment equation for Hookean dumbbells in a ^-solvent. I It 
has the analytical solution,! 

(Q*Q*) e = 1 - f* ds* e-C*-**) 7[0] (**, s*) (56) 

J — oo 

where, 7[ ] is the codeformational relative strain tensor.i The first order second 
moment equation has the form, 

^(Q*Q*) Z = k*-(Q*Q*) z +(Q*Q*) z -k* t -(Q*Q*) z +^ (e~ (Qt2/ ^ 2) Q*Q*) e (57) 

The ^-solvent distribution function ipg is a Gaussian,! and consequently, the complex 
moment on the right hand side of equation (|57| ) can be reduced to a function of 
(Q*Q*) e . The following equation is obtained on performing this reduction, 

■^(Q*Q*) Z = «• ■ (Q*Q*) Z + (Q*Q*) Z ■ ** T - (Q*Q*) Z + y (58) 

where, 

Y = , r/ * l(Q*Q*) e + 1 ] _1 ' (Q*Qle 

y/det [ (Q*Q*)e + /" 2 1] 

Clearly, equation fl58p could have also been derived by expanding the second 
moment equation for the Gaussian approximation fl42|) to first order in z. It follows, 
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therefore, that the Gaussian approximation is exact to first order in z. This is also 
the situation in the case of the Gaussian approximation introduced for the treatment 
of hydrodynamic interaction effects, where it was found to be exact to first order in 
the strength of hydrodynamic interaction, 

Equation (^) is a system of linear inhomogeneous ordinary differential equations, 
whose solution is, 

(Q*Q*) z = [* ds*e- {t *- s * ) E{t*,s*)-Y -E T {t*,s*) (59) 

where, E is the displacement gradient tensor.i 

The expression for the stress tensor (0) can also be expanded to first order in z. 
After reduction of complex moments to second moments, the stress tensor depends 
only on the second moments (Q*Q*) e and (Q*Q*) Z - Equations ( |56"1) and (|59|) may 
then be used to derive the following first order perturbation theory expression for 
the stress tensor in arbitrary homogeneous flows, 

-£- f = Y + f^ds*e-^~ s ^ (j m (t*,s*)-E(t*,s*)-Y.E T (t\s*)) (60) 

Note that Y, the direct contribution to the stress tensor, is isotropic only in the 
limit fi — > 0. 

The form in steady shear flow, of the tensors "yr i and E, has been tabulated in 
reference. @ Using the expression for the stress tensor (|60"1), and the definition of the 
viscometric functions (|TTD, the following first order perturbation theory results for 
the viscometric functions are obtained, 

* - 1 + (L^m x (61) 



X H nk B T \ v/T+T^A 3 / 2 

X 2 H nk B T + I v/rT7?A3/2 



(62) 



where, A = (1 + /i 2 )[l + /i 2 + 2A^7 2 ] — Xjj'j 2 - In the limit of A#7 going to zero, 
the expression for the viscosity fl6~l|) reduces to the expression derived earlier in 
section 5.1 — using the retarded motion expansion and the equilibrium perturbation 
expansion — for the zero shear rate viscosity. One can also show that tr(Q*Q*) 
reduces to the equilibrium moment q± [see equation fl38|)], in the limit Xui — > 0. 

The first order perturbation results are compared with simulation results, and 
with results of the Guassian approximation, in section 7. 



6 Fixman's Theory 

Many years ago, in path-breaking seminal work, Fixmani considered the simulta- 
neous inclusion of hydrodynamic interaction and excluded volume in bead-spring 
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chain models for dilute polymer solutions. In order to render the problem solvable, 
Fixman introduced a number of approximations. Since we are only concerned with 
excluded volume in the context of Hookean dumbbells in this work, we shall only 
consider those approximations which are relevant in this context. The introduc- 
tion of the quadratic potential, the governing equations of Fixman's theory, and the 
calculation of material functions predicted by the theory in simple shear flow are 
considered in this section. 



6.1 The quadratic potential 

With regard to excluded volume, the most crucial approximation of Fixman! is the 
replacement of the delta function potential with a quadratic potential. By adopting 
a Boson operator formulation of the governing equations, Fixman has shown that 
the delta function potential ( p7|) may be represented by, 

E*(Q*) = ^Q* G* Q* (63) 

where G* is a symmetric function of various configuration dependent quantities 
introduced in the Boson operator formalism. From this expression it is clear that 
a quadratic potential for the excluded volume may be obtained by replacing the 
fluctuating quantity G* with an average. Fixman obtains a quadratic potential by 
replacing G* with a configuration dependent average, as described below. 

As a result of replacing G* with its average, one can show from equation fl6"3"| ) 
that, 

f) f)E* 

^ = ! ~WW ) (64) 

In other words, for any given potential E*, one can find (G*) provided that the 
non-equilibrium distribution function ip(Q,t) with which to carry out the average 
on the right hand side of equation (|64[) is known. It turns out that ip(Q,t), in the 
presence of a quadratic excluded volume potential and with consistently averaged 
hydrodynamic interactionjlllll is a Gaussian distribution. Evaluation of the average 
for a 5-potential (|27|), with a Gaussian distribution (f|l]), leads to the following 
non-dimensional quadratic potential, 

E*(Q*) = -- ; Z Q* (Q*Q*y 1 Q* (65) 
2 Jdet(Q*Q*) 



It is straightforward to show that the potential ( pop leads to an unphysical non- 
central excluded volume force between the beads. Fixman, perhaps for this reason, 
introduces a further approximation which consists of replacing the above potential 
with the following simpler form, 

1 z O* 2 

E*{Q*) = --- g (66) 

2 « 2 yJdet(Q*Q*) 
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where, a is defined as before by equation (f45[) . However, in this case, a obeys the 
consistency relation, 

z = (a 2 -l)a 3 (67) 



Interestingly enough, equation (|6|) reduces to equation (l67f) in the limit // — > 
(which corresponds to a (^-potential). 

It is appropriate here to note that while Fixman has presented all his arguments 
for bead-spring chain models, the form of Fixman's potential for dumbbells given 
above can be found in the book by Larson.0 

The consequences of adopting the quadratic excluded volume potential ( |66"1) are 
briefly discussed in the following section. It is worthwhile to point out here that 
Fixman's original formulation of the problem was not in terms of the Gaussian 
distributions and second moments discussed below. Rather, his attempt was to 
directly solve the diffusion equation (P for the configurational distribution function, 
and then carry out the average in equation (g) to obtain the rheological properties. 
It is, however, possible to discuss his approach within the framework developed 
subsequently by Ottinge and this is the procedure that is adopted here. A graphic 
exposition of Fixman's algorithm is given in reference.0 



6.2 The governing equations 



With excluded volume interactions described by the quadratic potential (|66"D, the 
diffusion equation (p]) becomes linear in the bead-connector vector. As a result, the 
diffusion equation is exactly satisfied by a Gaussian distribution pT|). In Fixman's 
theory therefore, a tractable model is obtained not by approximating the distribution 
function, as in the case of the Gaussian approximation, but by introducing the 
quadratic potential (|66|). 

While ip{Q, t) is a Gaussian both in the Gaussian approximation and in Fixman's 
theory, the second moment which completely determines these distributions is dif- 
ferent in the two cases. In Fixman's theory, the non-dimensional second moment 
(Q*Q*) is governed by the equation, 



-^(Q*Q*) = k*-(Q*Q*)+(Q*Q*)-k* 



l - 



a 2 Jdet{Q*Q*} 



(Q*Q*) + 1 (68) 



Note that the second moment equation reduces to equation (|67|) at equilibrium. 
The Giesekus expression for the stress tensor has the form, 



nk B T 



a 2 Jdet(Q*Q*) 



(Q*Q*) + 1 



(69) 



As a result, the stress tensor in Fixman's theory, for any flow situation, may be 
obtained once equation (|I8|) is solved for (Q*Q*). 
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It is clear from equation fl69|) that in Fixman's theory, as for /x > in the Gaussian 
approximation, rheological properties are affected both directly and indirectly by the 
presence of excluded volume. 

As will be shown in the section below, at steady state in simple shear flow, the 
problem of solving the governing equation for the second moments fl68|), reduces to 
one of solving a single nonlinear algebraic equation. In the linear viscoelastic limit, 
however, analytical expressions for the various properties can be derived in the same 
manner as described earlier for the Gaussian approximation. Indeed, it can be shown 



that the linear viscoelastic properties are given by equations flag ) and (54), where 
the quantities Tt and r* are now given by, 

H = 1 ; t* = a 2 (70) 

In steady simple shear flow, substituting equation Q55D for (Q*Q*) and equa- 
tion flTOT) for k, into the second moment equation (|68|), leads to the following equa- 
tions for the components of (Q*Q*), 

si = (1 + 2 \ 2 H j 2 si ) s 2 ; s 3 = s 2 ; s 4 = \ H A fsl 

where, s 2 must satisfy the nonlinear algebraic equation, 

sT - sT = ^ = 1} (71) 

The normalised viscometric functions can be found by using the definitions (|TTj) , 
and the results above for the zero shear rate properties, 

77 Pi o a 2 ^i,o « 4 

The nonlinear algebraic equation ([71]) is solved here with a Newton-Raphson 
scheme. The material functions predicted by Fixman's theory are compared with 
the predictions of the narrow Gaussian potential in the section below. 



7 Results and Discussion 

The predictions of rheological properties in simple shear flow, by the various theories 
for the excluded volume effect, are compared in this section. Predictions in the limit 
of zero shear rate are first considered below, and those at finite non-zero shear rates 
subsequently. 

7.1 Zero shear rate properties 

Figure [I] is a plot of (%,,o/ uksT) versus /i for z = 3 and z = 100. The con- 
tinuous curves are exact predictions obtained by numerical quadrature, while the 
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Figure 1: Non-dimensional zero shear rate viscosity versus the extent of excluded volume 
interaction /i, for two values of the strength of the interaction z. The continuous lines are exact 
predictions obtained by numerical quadrature, the triangles and circles are results of Brownian 
dynamics simulations, the dashed and the dot-dashed lines are the approximate predictions of the 
Gaussian approximation, and the first order perturbation theory, respectively, and the filled squares 
are the predictions of Fixman's theory. The error bars in the Brownian dynamics simulations 
cannot be resolved within the line thickness. 

triangles and circles are exact results of Brownian dynamics simulations carried out 
at equilibrium (without variance reduction) [see equation ([Ti])]. The dashed lines 
are the predictions of the Gaussian approximation for the narrow Gaussian potential 
[i.e. equation (0), with r* and TC given by equations ([||) and (|50"D , respectively], 
the dot-dashed curve is the prediction of the first order perturbation theory [i.e. 
equation (|61|) in the limit Xnj — ► 0], and the filled squares are the predictions of 
Fixman's theory [i.e. equation (0), with r* and H given by equations ([70])]. The 
parameter /i does not enter into Fixman's theory, however, these values are plotted 
corresponding to \i — 0, since the quadratic potential is used in Fixman's theory as 
an approximation for the (^-potential. 

The first feature to be noticed in figure |l] is the reassuring closeness of the exact 
results obtained by using the retarded motion expansion and by Brownian dynamics 
simulations. The retarded motion expansion provides a means of validating the 
results of Brownian dynamics simulations. 

In the limit \i — *■ 0, and for large values of yU, the continuous curves and Brownian 
dynamics simulations reveal that, as expected, the exact predictions of the narrow 
Gaussian potential tend to the z-independent ^-solvent value, {r] p fl/\ H nk B T) = 1. 
This implies, as pointed out earlier, that the use of a 5-function potential to represent 
excluded volume interactions does not lead to any change in the zero shear rate 
viscosity prediction. On the other hand, figure |l] seems to suggest that a finite 
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range of excluded volume interaction is required to cause a change from the 6- 
solvent value. Away from these limits, at non-zero values of //, the narrow Gaussian 
potential predicts an increase in the value of zero shear rate viscosity. The existence 
of shear thinning in good solvents can be attributed to this increase. This follows 
from the fact that at high shear rates, as the effect of the excluded volume interaction 
diminishes, the viscosity is expected to return to its ^-solvent value. We shall see 
later that this expectation is indeed justified. 

The dashed lines in figure [l] indicate that in the limit of zero shear rate, for a 
given value of z, the Gaussian approximation is reasonably accurate above a certain 
value of /i. This limiting value of fi appears to be smaller for smaller values of z. 
A similar behavior is also observed with regard to the prediction of the zero shear 
rate first normal stress difference (see figure 0). Thus it appears that the exact 
configurational distribution function if>(Q, t) becomes increasingly non-Gaussian as 
the narrow Gaussian potential becomes narrower, and as the strength of the excluded 
volume interaction becomes larger. 

For the large values of z considered in figure |], results obtained with the first or- 
der perturbation expansion in z cannot be expected to be accurate. It is clear from 
the dot-dashed lines that the perturbation expansion results deviate significantly 
from exact results for small values of //. However, they become increasingly accu- 
rate as \i increases, for a given value of z. This can be understood by considering 
equation (|6lD, which indicates that in the limit A#7 — > 0, the first order correction 
to the 6*-solvent value increases as z increases, but decreases as \x increases. 

As n increases from zero, values of the zero shear rate viscosity predicted by the 
Gaussian approximation, approach the exact values more rapidly than the predic- 
tions of the first order perturbation theory, for a given value of z. The Gaussian 
approximation is a non-perturbative approximation — however, it was shown in sec- 
tion 5.4, to be exact to first order in z. One way to understand this is to consider 
the Gaussian approximation to consist of an infinite number of higher order terms, 
whose nature is unknown. In this sense, it is an uncontrolled approximation, which 
remains accurate at values of z and where the first order perturbation expansion 
becomes inaccurate. As will be seen shortly, these remarks apply also to the results 
obtained at finite shear rate. 

The difference in the prediction of the zero shear rate viscosity by Fixman's 
theory and by the narrow Gaussian potential is evident in figure [I]. It is also worth 
noting that, although the relaxation weight and the relaxation time are different in 
Fixman's theory and in the Gaussian approximation for \i — 0, they lead to the 
same prediction of the zero shear rate viscosity. This is, however, not true for the 
zero shear rate first normal stress difference. The ratio Uq v , defined by,li3 

^ = (73) 
% 

is equal, in the zero shear rate limit, to (2/a 2 ) in the Gaussian approximation, while 
it has a constant value of 2 in Fixman's theory. 
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Figure 2: Non-dimensional zero shear rate first normal stress difference versus //, for two values of 
z. The continuous lines are exact predictions obtained by numerical quadrature, the triangles and 
circles are results of Brownian dynamics simulations, the dashed lines are approximate predictions 
using the Gaussian approximation, and the dotted lines are approximate predictions using the 
uniform expansion model. The error bars in the Brownian dynamics simulations cannot be resolved 
within the line thickness. 

Both the uniform expansion model and the Gaussian approximation use Gaussian 
distributions in order to evaluate averages. However, the uniform expansion model 
uses different Gaussian distributions for different equilibrium averages, such that the 
best approximation is obtained. While this does not lead to any difference in the 
prediction of the zero shear rate viscosity by the two approximations, figure |2| reveals 
that, at small enough values of /x, there is a significant difference in the prediction 
of the zero shear rate first normal stress difference. Clearly, the uniform expansion 
model continues to be a reasonable approximation for values of fi at which the 
Gaussian approximation is no longer accurate. However, even the uniform expansion 
model leads to a poor approximation at sufficiently small values of [i. 

In the discussion of the equilibrium perturbation expansion in section 5.1.2, it 
was pointed out that, in principle, the equilibrium moment q\ can be obtained for 
any non-zero value of fi, provided that enough numbers of terms in the series for 
the quantities So and Si, are summed [see equations (|37D]. Figure 3 displays the 
sum So, for different numbers of summed terms in the perturbation expansion, as a 
function of /i, for z = 3.0. Exact results are obtained by noting that So = I , 
and evaluating J by numerical quadrature. Clearly, more terms of the expansion 
are required for convergence as fi decreases. The terms of the series, which are 
alternating in sign, keep increasing rapidly in magnitude, until n is approximately 
> (z/fi 3 ), before they begin to decrease. Therefore, as z increases, or as \x decreases, 
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Figure 3: The sum Sq versus fi, for different numbers of summed terms k in the perturbation 
expansion [see equation (5?])]. The exact results are obtained by numerical quadrature. 

more and more terms are required for the sum to converge. Above a threshold 
value of {z/fi 3 ) however, it becomes impossible to evaluate Sq since the round-off 
errors due to the summation of large numbers makes the perturbation expansion 
numerically useless. A similar problem is also encountered while evaluating the sum 
Si. In short, the hope of extrapolating finite /i results to the limit /x = 0, cannot be 
realised. Therefore, in the case of Hookean dumbbells, one cannot obtain equilibrium 
moments for a delta function excluded volume potential by using a narrow Gaussian 
potential and considering a perturbation expansion in the limit /x — > 0. 

7.2 Steady state viscometric functions 

The results of Brownian dynamics simulations (without variance reduction) dis- 
played in figure f| reveal that the dependence of the viscosity and the first normal 
stress difference on /i, at a value of the non-dimensional shear rate A#7 = 0.3, is 
similar in shape to the dependence observed in the limit of zero shear rate. At small 
and large values of \x the material functions tend to the ^-solvent value, and exhibit 
a maximum at some value in between. Since, even at this non-zero value of shear 
rate, it appears that the viscosity and the first normal stress difference remain at 
their ^-solvent values for \x = 0, it implies that the use of a (^-potential to represent 
excluded volume interactions would not predict any shear thinning. On the other 
hand, as we shall see subsequently, a quadratic potential does predict substantial 
shear thinning. At Xnj = 0.3, the Gaussian approximation seems to be accurate 
above roughly the same values of /x as were observed in the limit of zero shear rate. 
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Figure 4: Non-dimensional viscosity and first normal stress difference versus \x for two values of 
z, at a non-dimensional shear rate A#7 = 0.3. The squares and circles are results of Brownian 
dynamics simulations, and the dashed and continuous lines are the predictions of the Gaussian 
approximation, for z = 3 and z = 100, respectively. The error bars in the Brownian dynamics 
simulations cannot be resolved within the line thickness. 

Figures || and |6] are plots of non-dimensional viscosity and first normal stress dif- 
ference versus the non-dimensional shear rate A#7. Figure [5] displays the dependence 
of these viscometric functions on the parameter /j,, for a fixed value of z — 0.1, while 
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figure H displays the dependence on the parameter z, for a fixed value of \i = 2.5. 
The prediction of shear thinning for non-zero values of /i is apparent, and our earlier 
expectation in this direction is justified. In particular, the predictions of Brownian 
dynamics simulations, the Gaussian approximation, and the first order perturbation 
theory tend to ^-solvent values at high shear rates. 

Shear thinning, which is seemingly physically meaningful, is also predicted, as 
can be seen from figure |^, by both the Gaussian approximation and the first order 
perturbation theory, for ji = 0. This corresponds to a 5-function excluded volume 
potential, and is clearly an artifact of the perturbation expansion, since rigorous 
calculations indicate a trivial result. It remains to be seen if the situation is different 
in the limit of long chains. 

For small enough values of z, and large enough values of /i, the results of the 
Gaussian approximation, and the first order perturbation expansion, agree exceed- 
ingly well with the exact results of Brownian dynamics simulations. Indeed, as A #7 
increases for fixed values of z and fi, both the Gaussian approximation, and the 
first order perturbation expansion become increasingly accurate. In particular, if 
both the approximations are accurate at zero shear rate, they continue to remain 
accurate at non-zero shear rates. This can be understood in the case of the first 
order perturbation expansion by considering equations (|6lD and flB^). Clearly, the 
departure from the ^-solvent values decreases as Ah7 increases. This is inline with 
the intuitive expectation of decreasing excluded volume interactions with increasing 
shear rate. 

For a fixed value of the shear rate A#7, as z increases, or \i decreases, the pre- 
dictions of the Gaussian approximation and the first order perturbation expansion 
become increasingly inaccurate, with the first order perturbation expansion breaking 
down before the Gaussian approximation. This can be expected, since the Gaussian 
approximation — being exact to first order in z — is at least as accurate as the first 
order perturbation expansion. 

The two different Brownian dynamics simulation algorithms mentioned in sec- 
tion 5.2 were used to obtain the data in figures |5] and || While the algorithm with 
variance reduction was used for shear rates upto A#7 = 0.1, the algorithm without 
variance reduction was used at higher shear rates. With regard to the results ob- 
tained by variance reduction, it was found that the variance was typically reduced by 
a factor of five to ten by the parallel equilibrium simulation subtraction procedure. 
While the magnitude of the reduced variance was relatively independent of shear 
rate for the viscosity, it decreased with increasing shear rate for the first normal 
stress difference. The time to reach steady-state, from start-up of flow, was roughly 
ten relaxation times for z = 0.1, z = 3 and z = 30, and roughly fifteen relaxation 
times for z = 100. Rheological properties in the two parallel simulations remained 
correlated during the time required to reach steady-state, and as a result the present 
technique proved adequate for the purpose of variance reduction. 

We have seen earlier — in figure [I] — that both Fixman's theory, and the Gaussian 
approximation for /i = 0, lead to identical values for the zero shear rate viscosity. 
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Figure 5: Non-dimensional viscosity and first normal stress difference versus non-dimensional 
shear rate A#7, for three values of fi at z = 0.1. The circles are results of Brownian dynamics 
simulations, the dashed lines are the predictions of the Gaussian approximation, and the triangles 
are the predictions of the first order perturbation expansion. In the viscosity plot, the error bars 
in the Brownian dynamics simulations are smaller than the size of the symbols. 



This coincidence is, however, restricted to the limit of zero shear rate. At non-zero 
shear rates, as can be seen from figure [F], there is considerable divergence between 
the predictions of the two theories. 
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Figure 6: Non-dimensional viscosity and first normal stress difference versus Xjjj for three 
values of z, at fi = 2.5. The symbols are as indicated in the caption to figure [5| The error bars in 
the Brownian dynamics simulations are smaller than the size of the symbols. 



Fixman's theory for dumbbells, though differing considerably from the narrow 
Gaussian potential in terms of its predictions of rheological properties, has the ap- 
pealing aspect that it captures some of the universal features of the behavior of 
good solvents — which can only be expected from bead-spring chain theories in the 
limit of a large number of beads. We have seen this universal behavior earlier 
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Figure 7: Non-dimensional viscosity versus non-dimensional shear rate for z — 100. The 
continuous line is the prediction of Fixman's theory, and the dashed line is the prediction of the 
Gaussian approximation for fi = 0. 

in the correct prediction of the end-to-end distance scaling, and in the parameter 
free nature of the ratio Uq, v . Figure |8] displays the prediction by Fixman's theory 
of the reduced variable (r) p /r) p> o) versus the non-dimensional shear rate (3 = A p 7, 
where, A p = ( [rj\o Mt] s /Na k& T ), is a characteristic relaxation time. Here, [77] o is 
the zero shear rate intrinsic viscosity, M is the molecular weight and N\ is Ava- 
gadro's number. For dilute solutions one can show that, (3 = r] Pi0 j /uk^T '. The 
figure clearly reveals that as z — > 00, the curves for different values of z overlap. 
Therefore, in this respect also, Fixman's theory mimics the universal behavior ex- 
pected of long chains. The source of this behavior can be understood by examining 
equation fl7iP . In the limit of a 1, one can show that S2 = [a 6 (Aj^) -2 ] 1 / 5 '. As 
a result, (r/ p /r/ P) o) = [a 2 A#7]~ 2 / 5 — leading to the observed scaling. At small values 
of A#7, (r} p /r} Pt Q) — > 1. One can therefore construct the analytical expression, 

HE- = (i + a 4\*rf)-W (74) 

and expect it to be accurate at very small and large values of (3. The dot-dashed 
curve in figure § shows that equation ([74]) is accurate over a fairly wide range of f3. 

The universal behavior discussed above is not exhibited by the Gaussian approx- 
imation. The reason for this can be easily understood in the limit /i — > 0, where, 
the second moment equation fl42|) reduces at steady-state to a non-linear algebraic 
equation. Indeed, one can show that for a ^> 1, (r/p/r/p.o) = [1 + A 2 ^ 2 ] -1 / 5 . As a 
result, the normalised material functions do not collapse onto a single curve when 
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Figure 8: Reduced viscosity versus reduced shear rate predicted by Fixman's theory for large 
values of z. The continuous line and circles are obtained numerically, while the dot-dashed line is 
obtained with the analytical expression (|74[). 



plotted versus (3. It is, however, not realistic to expect a dumbbell model to ex- 
hibit universal features, and the real verification of universal behavior requires the 
development of a theory for long bead-spring chains. 

The non-dimensional ratio U^ v [see equation fl73]) l has a constant value of two, 
independent of shear rate, for Hookean dumbbells in ^-solvents, and in Fixman's 
theory for good solvents. Figure 9 displays the predictions of by the narrow 
Gaussian potential, obtained by Brownian dynamics simulations, the Gaussian ap- 
proximation, and the following first order perturbation expansion (which can be 
derived from equations flST]) and (|62|) ), 



u ^ = 2 - 2z 7TtW^ (75) 



where, A has been defined below equation (p2|). Since a logarithmic scale has been 
chosen for the shear rate axis, it is difficult to represent the zero shear rate value 
of Uq/rj. However, since it is very nearly constant at low values of shear rate, the 
zero shear rate value is represented in figure 9 by the filled circles on the y-axis. All 
the data in figure 9 have the same trend of remaining nearly constant at low shear 
rates, and approaching asymptotically the value of two at high shear rates. In the 
case of the first order perturbation expansion, this can be understood by considering 
equation ( ff5| ) in the limit A#7 — > oo. 

The dependence of the mean squared end-to-end distance (Q* 2 ) on A#7 is re- 
vealed in Figure 10. The circles and triangles are the results of Brownian dynamics 
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Figure 9: The ratio [see equation (u3j)] versus The circles are results of Brownian 

dynamics simulations, the dashed lines are the predictions of the Gaussian approximation, and the 
dotted lines are the predictions of the first order perturbation expansion. The filled circles on the 
y-axis represent zero shear rate values of U^ n obtained by equilibrium simulations. The error bars 
in the Brownian dynamics simulations are smaller than the size of the symbols. 



simulations obtained with the algorithm without variance reduction. The dotted 
and dashed lines are plots of the following expression, 



(Q* Z )=3 + 2X 2 H Y + 



v/IT7?a 3 / 2 



3(1 + fi 2 ) + (5 + 6^) X 2 H Y + 2 A^7 4 (76) 



obtained from the first order perturbation expansion, and the continuous line is the 
well known result for a theta solvent,! (Q* 2 ) = 3 + 2\ 2 H A f 2 . Interestingly, the effect 
of excluded volume interactions on swelling increases with increasing shear rate. In 
the case of the first order perturbation expansion — which appears to be accurate 
for z — 3, but not for z = 30 — this can be understood by noting, in the expression 
above, that the term representing the correction to the theta solvent result due to 
the presence of excluded volume increases as the shear rate increases. 



8 Conclusions 

The use of a narrow Gaussian potential, to describe the excluded volume interactions 
between the beads of a Hookean dumbbell model, leads to the prediction of swelling 
and shear thinning for relatively small non-zero values of the extent of interaction 
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Figure 10: Non-dimensional mean squared end-to-end vector versus non-dimensional shear rate 
A#7, for two values of z at \i = 1. The circles and triangles are results of Brownian dynamics 
simulations, the dotted and dashed lines are the predictions of the first order perturbation expan- 
sion, and the continuous line is the analytical solution for a theta solvent. The error bars in the 
Brownian dynamics simulations are smaller than the size of the symbols. 



[i. This is essentially caused by an increase in the magnitude of the equilibrium mo- 
ments relative to their ^-solvent values. A delta function description of the excluded 
volume potential, on the other hand, is found to predict neither swelling nor shear 
thinning for Hookean dumbbells. 

For a given strength of the excluded volume interaction z, the Gaussian approxi- 
mation is found to be reasonably accurate for values of \x larger than some threshold 
value. The behavior of the Gaussian approximation can be understood by compar- 
ing its predictions with those of a first order perturbation expansion in z, since it is 
shown here to be exact to first order in z. The perturbation expansion reveals that 
the departure of the viscometric functions from their ^-solvent values increases with 
increasing z, but decreases with increasing //, and increasing shear rate XhI- 

The use of a quadratic potential in Fixman's theory leads to the prediction 
of viscometric functions which are considerably different from those of the narrow 
Gaussian potential. However, Fixman's theory for dumbbells reproduces a number 
of universal features observed in good solvents. 
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A The uniform expansion model 

The equilibrium average of any quantity X{Q) is given by, 

POeq = ^cq / X (Q) ^P{-^Q 2 ~ E *} d Q (77) 

where, b 2 = 3(fceT/if). By multiplying and dividing the integrand with the Gaus- 
sian distribution ?p' eq (Q) [defined by equation (|39|)1, and using the normalisation 
conditions on ip eq (Q) and ip' (Q), equation (|77| ) can be rewritten in the form, 

(X(Q) e- B( Qh' 



where, B(Q) = (3/2) [ (1/b 2 ) - (1/V 2 ) } Q 2 + E*. For ip' cq (Q) to be a good approx- 
imation to ip e<i {Q), B(Q) must be small. On expanding exp(— B(Q)) in a Taylor's 
series, one can then write equation ([78]), to first order in B(Q) as, 

(x) eq = (xy cq -(xBy cq + (xy cq (By eq ( 79 ) 

If V is chosen such that (X B)' eq = {X} eq (B) , eq , it then follows that (X) is well 
approximated by (X) . 
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The choice of b' that best optimises the approximation clearly depends on the 
quantity X(Q) that is averaged. In our case we require averages of the quantities 
X( m ) = Q 2rn , m = 1,2, 3. From the well known result for the moments of a Gaussian 
distribution, 



(Q 2m y 



9 r2K ri<n 



cq 



r(m + -) 



(80) 



where, b'^ represents the value of b' corresponding to Q 2m . It follows that the 
optimum value of b'^ is found by solving the following non-linear algebraic equation 

for a (m) , 



(a 2 



(m) 



1) m + 



{ 



2m 



Km) + ^] 3/2 + 



1} = 



(81) 



(m) 



where, a^ m ) = b', m Jb; m = 1,2,3. In terms of c*( m ), 



U, 



<) n(2p+!) 

P =i 



(82) 
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